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Abstract 

In this paper, we put forward to another route realizing topological su- 
perfluid (TS). In contrast to conventional method, spin-orbit coupling 
and external magnetic field are not requisite. Introducing an experimen- 
tally feasible technique called on-site rotation (OSR) into p-band hon- 
eycomb optical lattices for spinless Fermi gases and considering CDW 
and pairing on the same footing, we investigate the effects of OSR on 
superfluidity. The results suggest that when OSR is beyond a critical 
value, where CDW vanishes, the system transits from a normal super- 
fluid (NS) with zero TKNN number to TS labeled by a non-zero TKNN 
number. In addition, phase transitions between different TS are also 
possible. 



PACS number (s): 67.85. Lm, 03.65. Vf, 74.20.-z 



1 Introduction 

Topological superfluid (or superconductor) (TS) has a full pairing gap in the bulk and is 
labeled by a non-zero integer topological invariant [H[2]- From the famous bulk-boundary 
correspondence such a topological integer ensures the existence of gapless excitations on the 
boundary of the system, in other words Majorana fermions (MF) [3] in vortex core of pairing 
order parameter. Roughly, MFs are neither fermions nor bosons but non-Abelian anyons 
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[I] and play an important role for the realization of fault-tolerant topological quantum 
computation (TQC) [5]. The application prospect of MFs makes TS become one of the 
hottest frontiers. 

In the condensed matter physics some practical two-dimensional systems have been 
theoretically proposed to realize TS [SI El El El QUI [HI 121 CE3] • In terms of these systems 
the entrance into TS requires subtle adjustment of Hamiltonian and it is very difficult in 
condensed matter physics, although MFs have been detected in InSb nanowires contacted 
with one normal (Au) and one superconducting electrode (NbTiN) [14]. In the light of 
the disadvantage for condensed matter, TS has been also suggested in cold Fermi gases 
owing to their many controllable advantages and operabilities. Following the successful 
observation p-wave Feshbach resonance (FR), Gurarie et al. |15| show that degenerate 
Fermi gases near a p-wave FR naturally give a concrete realization of TS. Zhang et al. 
|16j propose to create TS directly from an s-wave interaction making use of an artificially 
generated spin-orbit coupling (SOC). In fact, SOC have been realized in a neutral atomic 
Bose-Einstein condensate (BEC) by dressing two atomic spin states with a pair of lasers 
and the same technique is also feasible for cold Fermi gases [17} [T8] . Realizing that in a 
dual transformation SOC is formally equivalent to a p-wave superfluid gap, Sato et al. [19] 
suggest to artificially generate the vortices of SOC by using lasers carrying orbital angular 
momentum. In terms of the latter two ways, SOC and a large magnetic field are crucial in 
order to enter into TS. 

In this paper we suggest to create TS from spinless Fermi gases in p-band honeycomb 
optical lattices with so-called on-site rotation (OSR), that rotates every lattice site around 
its own center but keeps the whole lattice intact and has been realized for triangular optical 
lattices |20j . As a matter of fact p-band Fermi gases in honeycomb optical lattices in 
absence of OSR have shown many interesting characteristics, such as ferromagnetism |21| 
and Wigner crystallization |22[ [23] associating with flat bands, f-wave superfluidity with 
conventional pairing interaction [21]. The motivation for this paper comes from the Wu's 
work on quantum anomalous Hall effect in the same system [25]. Under single particle 
picture, Wu found that an arbitrary non-zero OSR not only breaks time-reversal symmetry 
and changes the topological properties of the system, but also drives a topological phase 
transition when OSR is beyond a critical value. Here we add on-site attraction interaction 
between p-band Fermi atoms into Hamiltonian and ask a question whether OSR can drive a 
phase transition into TS. The results are positive and OSR brings phase transitions not only 
from normal superfluid (NS) to TS, but also among different TS. From another perspective 
our work also can be considered as an extension to [25j , where f-wave superfluidity without 
OSR is discussed. Thus we also investigate the effects of OSR on f-wave superfluidity. 

Experimentally the route to realize TS suggested here is also feasible. On the one 
hand by placing two electro-optic modulators at two of three laser beams which coherently 
superpose to form a honeycomb lattice, OSR is available as illustrated in |26j . On the other 
hand due to Pauli exclusion principle the occupation of p-band is very convenient as long 
as the lowest s-band is fulfilled. In addition, on-site attraction interaction can be enhanced 
by using atoms with large magnetic moments, such as 167 Er with m = 7fiB on which laser 
cooling has been performed [27j . In contrast to |16[ 119] . where a pair of extra lasers and 
a large magnetic field are needed to produce an effective SOC and split two SOC bands 
respectively, our system is much simpler. 



2 



The organization of this paper is as follows. In section 2, we give the model and at 
the mean-field level investigate the ground state of the system by numerically minimizing 
the thermodynamic potential. In section 3 by calculating TKNN number Itknn [25] of 
occupied bands addressing the topological properties of the model, the topological phase 
diagram is obtained. In addition we also investigate the properties of edge states to prove 
our results. A brief conclusion is given in section 4. 



2 Model and Mean-Field Ground State 

The honeycomb optical lattice was realized experimentally by using three laser beams with 
co-planar propagating wavevectors quite some time age [29]. It is well known that a hon- 
eycomb lattice is not a Bravais lattice and there are two inequivalent sites in a unit cell, 
denoted by A and B respectively. Fulfilling the lowest s-band and defining three unit vectors 
e\ = n^-e x + \&y, e*2 = — &x + 5% and £3 = —e y , the Hamiltonian of p-band honeycomb 
optical lattices with OSR is 

H = t \\ Y [p?,iPr+Si,i + H.C. -U Y p\ x P\yPryPrx ~ ^ Y ~ V Y ™* W 

reA,i ' reA(BB r<=A<SB r£A®B 

where pp^ = (p^ x e x + p^ y e y ) ■ e*j and p? )X (p^ y ) is the annihilation operator for p x (p y ) 

band at the lattice site r. h ? = p\ x p?, x +P^ y Pr, y and lp >z = -i{p^ x Pr,y ~ p\ >y Pr,x) represent 
particle number and orbital angular moment operators. t\\ is the nearest-neighbor hopping 
matrix element of atoms in a bonds and positive due to the odd parity of the p-orbital. U 
(> 0), O (> 0) and [i are the on-site interaction strength, on-site rotation angular velocity 
and chemical potential, respectively. Note that we have neglected the nearest-neighbor 
atom hopping of it bonds and supposed the nearest neighbor distance in the lattice to be 
unit. 

When U = 0, introducing operator <p{k) = \pAx(k) , PAy{k) , PBx(k) , PBy{k)] T and making 
a unitary transformation <f) n (k) = U nm (k)^> m (k) , Hamiltonian can be diagonalized exactly. 
Meanwhile four energy bands can be obtained. Wu found two of four bands always are 
topological for any nonzero OSR and the others can be topological only if OSR is beyond 
a critical value [25] • On the basis of this findings, Wu proposed an orbital analogue of the 
quantum anomalous Hall effect, arising from orbital angular momentum polarization due 
to OSR. With O = 0, Lee et al. discussed f-wave superfluidity and charge density wave 
(CDW) in this system at the mean-field level |24j . Their results show that away from the 
half filling the system is f-wave superfluidity, while around the half filling superfluidity and 
CDW coexist and the system is a supersolid. Although superfluidity exists all the time, it 
is not topological as stated below. 

Following the same spirit in [24J we decouple interaction term into CDW channel 



Y (~u - ^)pi T P,, T + Y(-j u + 

rSA r£B 



tjCDW _ \^ 
J^int — / , 
T=x,y 

and pairing channel 

HfT n9 = - Y [^AP ] Ax (k)p\ y {-k) + A B pUk)p By (-k) + H.C. 



(2) 
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A nm {k')tf n [k')& m {-k') + F.C.I (3) 



where n =< hp A + h? B > /2 is filling factor of every site, Acdw = U < hf A — n? B > /2, 
Aa = UY,k < PA y (-k)pAx(k) >, A B = UY,k < PBy(-k)pBx(k) > are order parameters 
for CDW and superfluidity. In ([3]) we also express the pairing channel using quasiparticle 
ty(k). In this representation 

A nm {k) = A A [U? n (k)UZ m (-k)-U* n (k)U? m (-k)} 

+ A B [U* 3n (k)UZ m (-k) - UUk)U* 3m (-k)} . (4) 

After the mean-field approximation, the Hamiltonian ([I]) becomes a BdG Hamilto- 
nian H = [<f><(k), 4>{— k)]Hk[(j>{k), (ft {— k)] T and the properties of system are completely 
decided by the 8x8 matrix H^. Diagonalizing H^, we attain spectrum ei(k) and corre- 
spondingly eigenvectors (fi(k) (i = 1,2,- • -,8). Due to particle- hole symmetry inherent 
in this BdG Hamiltonian, the spectrum are symmetric about zero energy and we assume 
€l (k) = -e s (k) > 0, e 2 (k) = -e 7 (k) > 0, e 3 (k) = -e 6 (k) > 0, e A (k) = -e 5 (k) > 0. Then the 
thermodynamical potential at zero temperature is 

F = \ E h 4 ^ " e iW - £ 2(fc) - e 3 (*!) - e 4 (k)} + § |Aa| 2 + ^\A B \ 2 + ^jA 2 CDW , (5) 

k 

where iV is the number of the unit cell. Below we numerically minimize thermodynamic 
potential F about A^, and Acdw for fixed interaction strength U. Without loss of 
generality we choose A^ to be real, As = \ As\e and U/tn = 3.0. 

Fig.l shows the solutions of the Hamiltonian ([T]) at the mean-field level for changing 
chemical potential fi and OSR £1. Due to the particle-hole symmetry we only concentrate 
on negative chemical potential. Fig. 1(a) describes the variation of Acdw- F° r ^ = 
CDW is robust, but when Q is beyond a critical value fi c it vanishes suddenly. This is 
due to the fact that the appearance of OSR changes the band structures of single particle 
and breaks the nesting condition for CDW. Numerically we find £l c /t\\ ~ 0.4 ~ 0.6 and is 
monotonically increasing as the function of chemical potential. Fig. 1(b) shows the effect 
of OSR on particle density and further exemplifies that the variations of band structures 
driven by OSR cause nonmonotonic behavior of particle density. In contrast, superfluid 
order parameters A^, A# are more interesting and shown in (c) and (d). On the one hand 
for O > f2 c , A^ = \Ag\ and with the increase of OSR superfluid order smoothly decreases 
until disappearance. This suppression mechanism of superfluidity consists in time-reversal 
symmetry broken caused by OSR. While on the other hand for Q < £l c Aa is still decreasing 
but As is increasing with O. In fact the increase of A^ originates from the redistribution of 
particle density between sites A and B, in other words the decrease of Acdw as seen in (a). 
Thus at the mean-field level our calculation suggests that (1) OSR weakens stabilities of 
CDW and superfluidity and (2) for Q < O c , superfluidity and CDW coexist and the system 
is a supersolid. 

The optimization of 6 leads to 9 = ir for all parameters we choose. Below we discuss 
the effects of OSR on pairing symmetry for > Q c . From [24J without OSR and away 
from the half filling the intraband pairings in (jl]) have f-wave symmetry with three nodal 
lines of k x = 0, k y = ±k x /\/3 and 7r/3 rotation symmetry [Fig.2(d)]. On introducing OSR, 
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in terms of the pairing magnitude, nodal lines degenerate into some disconnected regions 
where intraband gap disappears, and tt/3 rotation symmetry retains. However real and 
imaginary parts of pairing break tt/3 into ir rotation symmetry. In Fig. 2(a) (b) and (c) as 
an example we show the magnitude, real and imaginary parts of An. 

3 Topological Phase Diagram and Majorana Fermion Modes 

In this section we discuss topological properties of Hamiltonian ([1]) . In terms of our system, 
it explicitly breaks the time-reversal symmetry due to OSR. Thus TKNN number Itknn 
plays a central role in deciding topological nature of the system [2S]. TKNN number is 
defined, by eigenvectors <fi(k) (i = 5, 6, 7, 8) corresponding to negative energy spectrm of the 
matrix Hf., into Itknn = -^p L J cPkTr dA, where A is a matrix one-form Ay = A^(k)dk u 

with Aij(k) = ip\(k)Vk v fj(k). By numerically calculating TKNN number [30], we show 
the topological phase diagram of the system in Fig. 3. For parameter region we choose, 
there are four different subregions labeled by Itknn = 1)0,-1,2 respectively. Moreover 
by comparison with Fig. 1(a) it is easily found that the boundary between Itknn = and 
other TKNN numbers in the direction of Q coincides with that of CDW disappearance. This 
finding is very important and ensures that topological order of our system is not topological 
CDW [31]. According to the criteria for TS [31 j Itknn = 2 corresponds to Abelian TS 
while Itknn = 1> —1 are non-Abelian TS. Thus Fig.3 tells us that OSR drives topological 
phase transition not only from NS to TS, but also between different TS. Here we mention a 
fact the energy gap of the bulk spectrum closes when topological phase transitions between 
topologically distinct phases occur. 

From the bulk-edge correspondence, a non-trivial bulk topological number implies the 
existence of gapless edge states localized on open edges of the system. Cold Fermi gases 
with sharp edges may be realized along the lines proposed in [32]. In order to understand 
the relation between Itknn and the number of edge states, we study the Hamiltonian ([TJ 
with the open boundary condition along the zigzag edge of the honeycomb lattice. The 
resulting excitation spectrum are depicted in Fig. 4 for representative parameter choices. 
Very explicitly the number of gapless states for every edge is one-to-one correspondence 
with the TKNN number. For Itknn = ±1 (Itknn = 2) there are one (two) pair(s) of 
gapless states, while for Itknn = 0, gapless state does not exist. Due to particle-hole 
symmetry, in terms of gapless states, they are Majorana fermion modes. It should also 
be remembered that the core of a vortex is topologically equivalent to an edge which has 
been closed on itself. The edge modes we describe are therefore equivalent to the Majorana 
fermions known to exist in the core of vortices of p-wave superfluids [33J . 

4 Conclusions 

In conclusion at the mean- field level we have investigated the effects of OSR on CDW and 
superfluidity for p-band spinless Fermi gases in honeycomb optical lattices. We found that 
OSR weakens the stabilities of CDW and superfluidity simultaneously, although superfluid- 
ity can survives a larger OSR. This conclusion leads to another important result that once 
CDW drops out the system enters into topological superfluidity. By numerically calculating 
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the TKNN number we obtained topological phase diagram of the system. In addition edge 
states, i.e. bulk-boundary correspondence are also investigated. 
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Figure 1: The mean-field solution of the Hamiltonian ([I]). Parameter U/tu = 3.0. 
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Figure 2: The symmetry of intraband pairing An. In (a) the magnitude, (b) real part 
and (c) imaginary part of An for fi/tii = 1.0 are shown. For comparison (d) plots An for 
fi/i|| = 0. 
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Figure 3: Topological phase diagram of the Hamiltonian ((TJ) at the mean-field level. The 
light grey, dark grey, black and blue colors correspond to Itknn = 1? 0, —1, 2 respectively. 
Parameter U/t\\ = 3.0. 
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Figure 4: The gapless edge states with the open boundary condition along the zigzag 
edge of the honeycomb lattice. In (a) Itknn = 0, fi/tn = —0.75, fi/iii = 0.3, A^All = 
0.969, AbA|| = 0.063, A CD w/t\\ = 1-462, (b) Itknn = 2, n/t\\ = -0.5, Q,/t\\ = 0.8, 
A A /t|| = AbA|| = 0.288, Acdw/H = 0, (c) I TKNN = -1, n/t\ { = -0.85, Sl/t\ { = 0.8, 
A A A|| = AbA|| = 0.523, A C dw/% = 0, (d) Itknn = 1, /"A|| = "0-65, O/t,, = 1.1, 
A A A|| = AbA|| = 0.307, AoowAll = 0. 
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